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ABSTRACT 

New data for calculation of the ionization and recombination rates have have been published in the 
past few years. Most of these are included in CHIANTI database. We used these data to calculate 
coUisional ionization and recombination rates for the non-Maxwellian ^-distributions with an enhanced 
number of particles in the high-energy tail, which have been detected in the solar transition region 
and the solar wind. Ionization equilibria for elements H to Zn arc derived. The K-distributions 
significantly influence both the ionization and recombination rates and widen the ion abundance 
peaks. In comparison with MaxwcUian distribution, the ion abundance peaks can also be shifted to 
lower or higher temperatures. The updated ionization equilibrium calculations result in large changes 
for several ions, notably Fe VIII-XIV. The results arc supplied in electronic form compatible with the 
CHIANTI database. 

Keywords: Atomic data - Atomic processes - Radiation mechanisms: non-thermal - Sun: corona - 
Sun: UV radiation - Sun: X-rays, gamma rays 



1. INTRODUCTION 

One of the most widely used assumptions in the in- 
terpretation of astrophysical spectra is that the emit- 
ting system is in thermal equilibrium. This means that 
the distribution of particle energies is at least locally 
Maxwellian, and can be characterized by the Boltzmann- 
Gibbs statistics which has one parameter, the tempera- 
ture. 

The generalizat ion of t he B o lzman n-Gibbs statis- 
Tsall isl (fl988l [20091) results in 



tics proposed by 
distributions (e.g., 



Leubncr 2003, l2004allbl : ICollied[200l 



iLivadiotis fc McComasI l2009il , characterized by a pa- 
rameter n and exhibiting a near-Maxwellian core and 
a high-en ergy power- l aw ta il (Sect. [21). First pro- 
posed by IVasvliuiiaa (|1968f l. the k distributions are 
now used to fit the observations of a wide vari- 
ety of astrophysical environments, e.g. in-situ mea- 
surements of particle distri butions in planetary mag- 
netic environments fe.g.. [ Pier rard fc Lemairc' "'""'" 



199C 



Mauk ct al. 2004; S chippers et al.. 2008: Xia o et al. 20Q8|; 
Dialvnas ct al. 200 3) and solar wi nd fe.g.. iCollier et alJ 



~996; .Maksimovic et all Il997al[bl : iPierrard et al I 



1999 : 



2009L 



Nieves-Chinchilla fc Vifiasl 120081: iLe Chat et alT 
20111 : iPierrardI l2012f ). as well as photon spec t ra of 
solar flare plasm as (e.g., iKasparova fc Karlickvl 120091 : 
lOka et aLl l2013f ). emission lin e spectra of plane tary 
neb ulae and galactic sources (jNichoUs et al.l2012[ : see 
alsolBinette et al.l2012 ) and even solar transition region 
(jDzifcakova et al.ll2011l ). A review on the K-distributions 



and their ap plications in astrophysic al plasma can be 
found, e.g., in IPierrard fc Lazair(|2010( ). 
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In the solar corona, presence of the K-distributions, 
or distributions exhibiting high-energy tails, can be ex- 
pected due to particle acceleration processes arising 
as a result of "nanoflare" heating. The nanoflares 
are an unknown energy release process of impulsive 
nature, occurin g possibly in storrns heating the so- 
lar corona (e.g.. iTripathi et al.ll2010t [Viall fc Klimchukl 
[20Tlt [Bradshaw et all 120121 : iWinebarged I2012D . While 
the direct evidence for enhanced suprat hermal popula- 
tions in the solar corona is still lacking ( Feldman et al.l 
[200l[Ha^ah et al.l[20Tol ). lPinfield eTalllligggl ) reported 
that the intensities of the Si III transition region lines 
observed by the SOHO/SUMER instrument do not cor- 
respo nd to a single Maxwellian dist ribution. Using their 
data, iDzifcakova fc Kulinoval (|2011| ) showed that the ob- 
served intensities can be explained by K-distributions 
once the photoexcitation is taken into account. These 
authors diagnosed k = 7 in the active region observed on 
the solar limb. Higher values of k were diagnosed for the 
quiet Sun and coronal hole, indiciating that the depar- 
tures from the Maxwellian distribution can be connected 
to the local magnetic activity. The diagnostic method 
also works for inhomogeneous plasmas characterized by 
differential emission measure. That the K-distributions 
can be present in the solar corona is also suggested 
by their presence in the so lar wind (e.g.. IPierrard et al.l 
[T999l : [Vocks fc Mannl[200l . 

Direct diagnostics of K-distributions in the solar 
corona using extreme-ultraviole t lines observed by 
the Hinodc/EIS s pectrometer (iCulhane et al.l 20071 ) 
were attempted bv IDzifcakova fc Kulinoval ()2010D and 
iMackoviak et al.l ()2013[ ). These authors proposed meth- 
ods for simultaneous diagnostics of the plasma tempera- 
ture, electron density, and k. However, majority of the 
line ratios sensitive to K-distributions suffer from poor 
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photon statistics, errors in atomic data and/or plasma 
inhomogeneities. In spite of this, one of the main resuhs 
of these works is that the sensitivity to ^-distributions, 
or to departures from the Maxwehian distribution in gen- 
eral, is enhanced if the line ratios involve lines originat- 
ing in neighbooring ionization stages. This is due the 
sensitivity of the line emissivity to the abundance of the 
emitting ion that depends directly on the ionization equi- 
librium in tu rn highly dependent on the type of the distri- 
bution fe.g..'D zifcakovall2002l : rWannawichian et al1l2003l : 
iDzifcakova 200^^ 

In the past decade, new calculations of the ioniza- 
tion, recombination rates and ionization equilibirium 
for the Maxwellian electron distribution were published. 
These are summarized in the continually updated CHI- 
ANTI dat abase, curre ntly available in version 7.1 (e.g., 
IDere et aLlHoOTl [20091: Itandi et al.l[20Tl [20T1 . These 
new calculations of the ionization and recombination 
rates result in significant differences with respect to the 
earlier calculations, e.g. of iMazzotta et al.l ()1998[ ). 

The availability of accurate atomic data are of cru- 
cial importance in correct determination of the proper- 
ties of the radiating astrophysical environment, and the 
solar corona in particular. In this paper, we present up- 
to-date calculations of the ionization and recombination 
rates (Sect. |31), and ionization equilibria (Sect. \^ for 
K-distributions for elements from H to Zn. Such calcula- 
tions are necessary for diagnostics of the K-distributions 
in both the solar transition region and the corona, as 
well as subsequent calculations of the radiative losses 
(jDudik et al.l 120 lit) or re sponses of variou s extreme- 
ultraviolet or X-ray filters (jDudik et al.ll2009l ) used both 
to model and observe these portions of the solar atmo- 
sphere, with potential applications to other astrophysical 
environments. 

2. THE NON-MAXWELLIAN k-DISTRIBUTIONS 

The K-distributions of electron kinetic energies £ rep- 
resent a family of non-Maxwellian distributions charac- 
terized by two parameters, k and T 



U£)d£ = A, 



£'/^d£ 



7ri/2(fcBT)3/2 (1 + 



(K-1.5)fcT 



)K+1 



(1) 



where the A = r(K+l)/ [r(K - 0.5)(k - 1.5)^/2] is the 
normalization constant, kB = l-i8 xl0~23 Jkg"-*^ is the 
Boltzmann constant, and kG (3/2,-|-oo), Tg(0,-|-c»). 
We note that the definition of K-distributions in Eq. ([T] 
top) corresponds to the ^-distributi ons of the second kind 
(e.g., iLivadiotis fc McComasll2009D . 

The shape of the distribution is controlled by the 
parameter k. Maxwellian distribution is recovered for 
K— >--|-oo, while the departures from the Maxwellian 
increase with k— >3/2. The departures from the 
Maxwellian distribution with decreasing k include in- 
crease of the number of particles in the high-energy 
tail as well as increase in the relative number of low- 
energy electrons (Fig. [T]fop). However, the mean energy 
{£) — 2)/2k-QT of the distribution does not depend on k. 
This allows for calculation of all quantities depending 
on the mean ene rgy of the distribution, e.g. pressure 
(|Dzifcakoval 120061 ) . The p arameter T has in the frame 
of nonextensive statistics (iTsallisI 119881 I2009D an anal- 
ogous meaning as thermodynamic temperature in the 



Boltzmann-Gib bs statistics. The reader is referred e.g. 
to the work of ILivadiotis fc McConiasI (|2009l I2010D for 
details. 

While the shape of the K-distribution differs from the 
Maxwellian with the same T at all energies E, the core of 
the K-distribution can be approximated by a Maxwellian 
distri bution with lower Tc = T{k — 3/2)/ k (|Oka et al.l 
|2013[ ). An example for k = 5 and T = 1MK is shown 
in Fig. [1] bottom, where the approximating Maxwellian 
distribution has Tq = 0.7T and has been scaled by the 
factor C 



C = 2.718. 



r(K + i) 



-3/2 



1 



1 



(2) 



r(K-i/2) 

(|Okaet al.l 120131 ) ■ which is «0.84 for k = 5. The differ- 
ence between the fniT) and /Maxw(T'c) are then mainly 
in the pronounced high-energy tail. This shows that 
the K distributions offer straightforward approximation 
of situations where a non-Maxwellian high-energy tail is 
present. 

3. IONIZATION AND RECOMBINATION RATES 

To calculate the ionization and recombination rates 
for K-distributions, we use the atomic data available 
through the CHIANTI database for a strophysica l spec- 
troscopy of opti cally thin plasmas (jDere et al.l 119971 : 
iLandi et al.ll2013l ). The analytical functional form of the 
^-distributions allows for relatively simple direct inte- 
gration of the ionization cross-sections fSect IX^ . How- 
ever, the recombination cross-sections are not contained 
in CHIANTI. These then have to be reverse-engineered 
from the Maxwellian recombination rates using assump- 
tions detailed in Sect. 13.31 

We note that the CHIANTI database allows for the 
treatment of non-Maxwellian distributions only if these 
can be represented by a series of individual Maxwellian 
distributions (with different Ts). The technique for cal- 
culation of ionization and recombination rates presented 
here can in principle be extended for any type of particle 
distribution, not only K-distributions. 

3.1. Atomic Data 

The CHIANTI atomic database since version 6 
(jDere et al.l |2009() contains continually updated ion- 
ization equilibrium for the Maxwellian distribution. 
This ionization equilibrium utilizes the cross-sections 
for direct ionization and autoionization, and_the corre- 
sponding rate coefficients from the work of lDerd ()2007D . 
Dielectronic and radiative recombination coefficients for 
the H to Al and Ar isoelectronic sequences are taken 
from the works of N. Badnell and colleagues, listed 

at http:/ / am.dpp . p hys. strath, ac. u k/ tam.oc /D A TA / 

(iBadnell et al.l (20031: IColgan et al.1 l2(M 12001 
Mitnik fc Badnelll l200i iBadnelll 120061: lAltun et all 



20051 I2OO6LI2OO 






iZatsarinnv et al.l 12005 
Bautista fc Badnelll 120071: 
Abdel-Nabv et al.l 120121 )" For the Si to Mn 
tronic sequences, the recombination rates 



iNikohc et aLT ' 



2006 



2010; 



isoelec- 

are based 

on the work s of iShull fc van Ste enberg (198 ^. | Nahai 



19961 I1997D iNahar fc Bauti sta' (2001|)r .Mazzotta et al 
19981) and iMazziteUi fc Ma ttioli (2002[K The reader 
3 referred to IDere et al.l ( 2009i ) for details. In the 



calculations presented in this paper, all atomic data for 
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ionization and recombination correspond to the ones 
used to produce the chianti.ioneq file in the CHIANTI 
v7.1. We note that ionization equilibria were pub hshed 
also by other authors (e.g. iBrvans et al.l l2009f ). but 
these use mostly the same atomic data as the CHIANTI 
database. 

3.2. Ionization Rates 

The ionization rates have been calculated by a numer- 
ical integration of the ionization cross section, ai, avail- 
able in CHIANTI database 



Ri 



:< aiV > 



2£ 
m 



1/2 



U£)d£, (3) 



for K = 2, 3, 5, 7, 10, 25, and 33. Each numerical integral 
is calculated as a sum of individual integrals splitted ac- 
cording to the ionization and auto-ionizaton energies. 

Typical changes in behaviour of the direct ionization 
and autoionization rates with K-distributions are shown 
in Fig. [2] for the ions C IV, Fe XII, and Fe XVII formed 
at transition region, quiet coronal, and flare conditions, 
respectively. For lower k, the temperature dependence of 
the ionization rates is much flatter, with lower maxima. 
Typically, for temperatures lower than those at which 
the ion abundance peaks (denoted by arrow in Fig. ^ 
see also Sect. HJ, the K-distributions result in increase of 
the ionization rates by up to several orders of magnitude 
with respect to the Maxwellian distribution. These de- 
viations from Maxwellian ionization rates increase with 
decreasing n. 

3.3. Recombination Rates 

Since the individual recombination cross-sections are 
not available in CHIANTI, the calculation of the re- 
combination rates for ^-distributi ons w as performed us- 
ing the method of iDzi fcakovai (|1992D , used also by 
iWannawichian et al.l ( 20031 ). This method allows for cal- 
culation of the recombination rates using approximations 
to the rates for the Maxwellian distribution. 

The cross-section ctrr for the radiative recombination 
is assumed to ha.v e a power-law dependence on energy 
(IOsterbrockl[T97l 



o'rr(^) ~ Ctxr/S 



»;+0.5 



(4) 



where Crr is a constant and 7/-I-0.5 is a power-law index. 
Subsequently, the radiative recombination rate for the k- 
distribution is of the form: 



OK 



4Crr r(K + 77-0.5)(K-1.5)T(1.5-?7) 



rr 



(27rm)i/2 r(K-0.5) (fcer)'' ' 

while the radiative recombition rate for the Maxwellian 
distribution is 



nMaxw ^ 4Crr r(1.5 - 1]) 

^^ (27rm)i/2 (fceT)') 

Therefore, it holds that 






R 



Maxw 
RR 



r(K + 7?-0.5)(K-1.5)'' 

r(K-o.5) 



(6) 



(7) 



The level of error introduced by this approximation is 
typically several per cent, i.e., lower than the error of 
the atomic data themselves. 



For the dielectronic r ecombination, foll owing approxi- 
mation has been taken ()Dzifcakovalll992f ) 



^DR — 



ri/3 ^ (1 + t,/{K - i.5)r)(«+i) 



(8) 



where parameters Oi and ti are the same as in similar 
expressions for the Maxwellian distribution 



R 



Maxw 
DR 



1 



7-1/3 Z^ 



V" Oi exp {-tjT) 



(9) 



where the coefficients a^ and ti arc provided within the 
CHIANTI database. The precision of this approximation 
is given by the magnitude of the second-order terms in 
the expa nsion of the coeffic ients a; and ti into series (Eqs. 
37-44 in lDzifcakoval[T99l . 

The typical behavior of the total recombination rates 
(^R,R, + ^dr) for K-distributions and for Maxwellian dis- 
tribution is shown in Fig. [21 It can be seen that the 
radiative recombination rate increases with decrease of 
K. This is a result of increasing number of low energy 
electrons for the ^-distributions (Fig. [1]), which domi- 
nate the recombination processes. The local change of 
slope of total recombination rate is caused by the con- 
tribution of dielectronic recombination. For some ions, 
e.g. Fe XVII, the contribution of dielectronic recombi- 
nation is dominant in temperature interval where these 
ions have a non-negligible abundance. 

4. THE IONIZATION EQUILIBRIUM 

Calculations of the collisional ionization equilibrium 
assume that there are no temporal variations in plasma 
temperature T. In the coronal conditions, the result- 
ing relative ion populations are given by the equilib- 
rium between the direct collisional ionization with auto- 
ionization and the radiative and dielectronic recombina- 
tion. Three-body processes can be neglect ed at low elec- 
tron densities typical in the solar corona (jPhillips et al.l 
I2008D . The radiative field is also usually assumed to 
be too weak, i.e., photoionization can be neglected as 
well. However, we note that photoionization may be im- 
portant for some transition-region ions with low ioniza- 
tion thresholds, but the effect will vary depending on the 
distance from the radiation field (e.g., the solar photo- 
sphere). 

The ionization equilibrium for K-distributions with 
K = 2,3, 5, 7, 10, 25, and 33 was calculated for ions of 
astrophysical interests, ranging fro m H (Z = 1 ) to Z n 
(Z = 30). Previo us cal culations of iDzifcakoyal (Il9 921. 
IDzifcakovai (|2002D . and IWannawichian et al.l (I200I in- 
volved only 12 most abundant elements. The current 
calculations are performed for the same set of tempera- 
tures as the calculations in the CHIANTI database in the 
chianti.ioneq file. I.e., the temperature spans the interval 
of log(r/K)e (4,9) with the step of Alog(T/K) = 0.05. 
The calculations are available in the same format as the 
chianti.ioneq file and are easily readable using the rou- 
tines read-ioneq.pro and plot_ioneq.pro available in CHI- 
ANTI running under SolarSoftware in IDL. More infor- 
mation on the CHIANTI .ioneq file format is provided 
in Appendix 1X1 

Examples of the current ionization equilibrium cal- 
culations for iron are in Fig. |3l and for carbon in 
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Fig. m The typical behavior is that the lower k, 
the flatter the ionization peaks. In addition, ioniza- 
tion peaks can be shifted to lower or higher T, de- 
pending on K, T and the individual ion. This means 
that a given relative ion abundance can be formed at 
a wider range of T for a K-distribution, with different 
peak formati on tempera t ure. S u ch behavior was al ready 
reported by iDzifcakoyal (I1992D . iDzifcakoval (|2002[) and 
iWannawichian et al.l ( 2003[ ). Typically, the shifts of Fe 
and C ion abundance peaks are to lower T with respect to 
the Maxwellian ionization equilibrium if log(r/K) ^ 5.5. 
At coronal temperatures, the Fe IX - Fe XVI ions are 
shifted to higher T. An interesting example is the ion 
Fe XVII, whose ionization peak is shifted to lower T for 
K lessapprox 3, but to higher T for k = 2 (Fig. ^bottom). 

There are a numbe r of differen c es wi th respect to the 
earlier calculations of IDzifcakoval (j2002[ ), which used the 
atomic data corresponding to those of IMazzotta et al.l 
()1998( ). These differences are illustrated in Fig. [S] In 
general, the lower the value of k, the greater the differ- 
ences with respect to the previous calculations. The most 
conspicuous examples are the Fe IX, which is shifted to 
slightly high er T instead to low er T as in the previous cal- 
culations of IDzifcakoval (|2002D . The Fe XII and Fe XIII 
ions are shifted to higher T (up to 0.1-0.15 dex). We note 
that the changes in the ionization equilibria due to the 
updated atomic dat a will reflect e.g. o n changes in the to- 
tal radiative losses (iDudik et al.ll201ll) and will also mod- 
ify the proposed dia gnostic m ethods for ^-distributi ons 
(jDzifcakovd fc Kuhnova 20iot JMackoviak et "alll2013[) . 

Ratios of relative abundances of individual ions can 
be used to diagnose the value of k and simultaneously 
T. S uch diagnostics can be a pplied e.g. in the solar 
wind (lOwocki fc Scuddeil I1983D . An example of diag- 
nostic diagrams is given in Fig. |6l We note that such 
diagrams cannot be directly applied on the observed line 
intensities unless the additional effect of K-distributions 
on the excitation and deexcitation rates is considered. 
However, these diagrams provide quantification of the 
expected changes to ratios of line intensities due to ion- 
ization equilibrium. Ratios of lines involving different 
ionization stages provide better options for determin- 
ing K, as noted already by (jDzifcakova fc Kulinoval[2010l : 
iMackoviak et al.ll2013h . 

5. SUMMARY 



Collisional ionization equilibrium calculations for re- 
distributions and optically thin plasmas were performed 
for all ions of elements H to Zn. To do that, the latest 
available atomic data for ionization and recombination 
were used. The calculations are available in the form 
of ionization equilibrium files compatible with the CHI- 
ANTI database, v7.1. 

For K-distributions, the ionization peaks are in general 
flatter and can be shifted to lower or higher T. This 
means that for K-distributions, individual ions are typ- 
ically formed in a wider range of temperatures, with 
different peak formation temperatures. Typically, ions 
formed at transition region temperatures are shifted to 
lower T, while the majority of coronal iron ions are 
shifted to higher T. Comparison to previous calculations 
is provided. Due to updated atomic data calculations, 
several ions in the present calculations are shifted to dif- 
fe rent temperatures with respect to previous calculations 
of lDzifcakoval ()2002D . The effect of k on the ratios of ion 
abundances is documented. 

The present calculations provide an accurate, neces- 
sary and useful tool for detection of astrophysical plas- 
mas out of thermal equilibrium, where the distribution 
of particles is characterized by an enhanced power-law 
tail. The supplied files compatible with the CHIANTI 
database and software should greatly facilitate syntheti- 
zation of both line and continuum spectra for optically 
thin astrophysical plasmas. 
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APPENDIX 
CHIANTI .lONEQ FILE FORMAT 

The Aoneq files, used by the CHIANTI database and software to store data on the ionization equilibrium, are in 
essence formatted ASCII files. An example of their format is given in Table [TJ The file begins with two numbers, 
N ^ giving the number of temperature points Tj, 1 < i < N , and .^max, denoting the maximum proton number Z for 
which the ionization equilibrium is provided. The next line gives the tabulated temperatures log(T!i/K). All other 
following lines begin with the Z and a positive integer -\-k denoting the roman numeral for the corresponding ion; e.g., 
26 and 9 stands for Fe IX. The line then contains relative ion abundances A^ '~^{T.i) in floating-point precision for the 
tabulated temperatures. 

Portion of the actual content of the kappa_05.ioneq file is given in Table [21 There, the relative abundances of the 
ions Fc IX-Fe XI arc fisted for k = 5 and for log(T/K) = 5.80 - 6.15. 
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Table 1 

ASCII format of the CHIANTI Aoneq file. Explanation of the symbols is given in the text. 







log(Ti/K) 


log(T2/K) 


log(r3/K) 


log(T4/K) 


log(T5/K) 


logCTej/K) . 


log(T;/K) 


los(T^f/K) 


1 
1 


1 
2 


A + '^iT,) 


At\T^) 


A + ^(T-3) 
aS^T3) 




A„ (T5) 


A + "(Te) . 
A + hT^) . 


■ a+\to . 


A + "(T„) 
A+l(Tjv) 


2 
2 

2 


1 
2 
3 




<(-2) 
Af^T,) 

■*H^(^2) 


AfiT,) 
A^tiTs) 


4 ("^' 

Ag(T4) 
A„;(T4) 


4J(T5) 


At>,, 

AtliTe) 

A+lln) ■ 


■ A+liTO ■ 


■ A + "(Tjv) 

■ A + 1(T„) 
A + 2(T„) 








Z 
Z 
Z 
Z 

z 


1 
+ (fc + 1) 

+ (Z + 1) 


A + "(Ti) 
A + ^Ti) 


A + ^Tj) 
A + ^in) 


Ap'iTs) 
A + ^Ta) 

-4 + ^(r3) 


Ap'lTi) 


a|"(T5) 
A+'=(T5) 


Af'iT,) . 

A+''{Te) 

A + ^{Ts) . 


. A + "(T., . 

A + ''{Ti) 
■ A + ^(T^) 


A + "(T„) 








Zmax 


+ (Z,„ax + 1) 
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Table 2 

Example of the kappa^OS.ioneq file produced for k = 5. The table shows relative abundances of the ions Fe IX - Fe XI. 



101 


30 






















5.80 


5.85 


5.90 


5.95 


6.00 


6.05 


6.10 


6.15 










26 
26 
26 
26 
26 


9 

10 
11 


0.357000 
0.0905900 
.. 0.00983700 


0.405000 
0.144700 
0.0228500 


0.422800 
0.211400 
0.0481900 


0.401500 
0.279200 
0.0912400 


0.342300 
0.329100 
0.153000 


0.257800 
0.340500 
0.223400 


0.167600 
0.302200 
0.277700 


0.091260 . 
0.223100 . 
0.284900 . 
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H TO Zn Ionization Equilibrium for the k-distributions 
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Figure 1. Top: The K-distributions with k = 2, 3, 5, 10, 25 and the Maxwellian distribution plotted for log(T/K) = 6.0. Bottom: The k = 5 
distribution and the approximation of its core with Maxwellian distribution with lower Tq. Colors and linestyles correspond to different 
values of n. (A color version of this figure is available in the online journal.) 
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Figure 2. Total ionization and recombination rates for C IV {top row), Fe XII {middle row) and Fe XVII {bottom row). Different linestyles 
correspond to different k. Black and red arrows denote maximum of the relative ion abundance for the Maxwellian and k = 2 distributions, 
respectively. (A color version of this figure is available in the online journal.) 
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Figure 3. Changes in the ionization cquihbrium with n for iron. Top: re = 10; Second row: k = 5; Third row: k = 5; Bottom row: k = 2. 
Individual ionization stages are indicated. (A color version of this figure is available in the online journal.) 
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Figure 4. Ionization equilibrium for Carbon. Different linestyles correspond to different k. (A color version of this figure is available in 
tfie online journal.) 
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Figure 6. Diagnostics of k from ratios of ion abundances. The values of k for each line are indicated. Thin lines represent isotherms for 
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